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Abstract 

We analyze the statistical properties of the spectrum of the QCD Dirac operator at low energy in a finite box of 
volume L 4 by means of partially quenched Chiral Perturbation Theory, a low-energy effective field theory based 
on the symmetries of QCD. We derive the two-point spectral correlation function from the discontinuity of the 
chiral susceptibility. For eigenvalues much smaller than m c = J ?2 /EL 2 , where F is the pion decay constant and 
E is the absolute value of the quark condensate, our result for the two-point correlation function coincides with 
the result previously obtained from chiral Random Matrix Theory (chRMT). The departure from the chRMT 
result above that scale is due to the contribution of the nonzero momentum modes. In terms of the variance of 
the number of eigenvalues in an interval containing n eigenvalues on average, it amounts to a crossover from a 
log n-behavior to a n 2 log n-behavior. 
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1 Introduction 



Since the seminal work by Wigner, Dyson and Mehta |l|, g, [3| on the Random Matrix Theory 
description of level correlations in nuclei, the problem of level statistics has been analyzed in great 
detail for many different quantum systems (see § for a recent review). In QCD, a statistical 
analysis can be applied to the eigenvalues of the Dirac operator ilp defined by 



= i\ k <p k . (1.1) 
The Euclidean QCD partition function for Nf quarks of mass rrif is given by 



, N f 

Z QCB = [ d A] u Y[det(il/) + m f ) e- SYM[A] , (1.2) 



and can thus be expressed as 

N f 

Z QCD = (JJJJ(iA fc +m/))YM. (1-3) 

f=l k 

Here, (• • -)ym denotes an average over the Yang-Mills partition function. This shows that, 
although the eigenvalues cannot be observed directly, their properties are of fundamental im- 
portance to the physics of QCD. 

For small enough energies, below the so-called Thouless energy [||, H fj], the eigenvalues 
of the QCD Dirac operator are strongly correlated, and their correlations are given by chiral 



Random Matrix Theory (chRMT) |j [K], 11]. For energy differences much larger than the 
Thouless energy but much smaller than Aqcd, the eigenvalues of the QCD Dirac operator 
show much weaker correlations that are different from chRMT. In this domain, the eigenvalue 
correlations can be computed perturbatively by means of partially quenched Chiral Perturbation 
Theory (pqChPT). This is a low-energy effective theory based only on the symmetries of QCD 
formulated to probe the spectrum of the QCD Dirac operator. Finally, for energies beyond 
-^■QCD the eigenvalues are uncorrelated. These different domains have been identified in lattice 



QCD simulations fl2| , 13, 14], and can be derived from random hopping models with the chiral 



symmetry of the QCD partition function p5| , |16| , 17, 19, [2(|. The occurrence of Random 
Matrix behavior in the Dirac spectrum can be understood naturally from universality arguments. 

The Thouless energy can be studied quantitatively by means of pqChPT. /,From the analysis 
of the average spectral density of the QCD Dirac operator, p(X), we found [@, |, |7| that the 
Thouless energy is given by m c = F 2 /T,L 2 , where F is the pion decay constant and £ is the 
absolute value of the quark condensate. The l/L 2 -dependence of m c is well-known from the 



theory of mesoscopic systems (see for example [28]). In addition to Aqcd there is one other 
important energy scale in the Dirac spectrum: the smallest nonzero eigenvalue A m j n . It is 
directly related to the spectral density near zero virtuality and is therefore given by the Banks- 



Casher formula [29]: A m i n = l/p(0) = Tt/Y>V . With the terminology adopted from the study of 



disordered mesoscopic systems |28j we can thus distinguish four different domains in the Dirac 
spectrum. The quantum domain and the ergodic domain are separated by A m i n , the ergodic 
domain and the diffusive domain are separated by the Thouless energy m c , and the diffusive 
domain and the ballistic domain are separated by Aqcd- 

In this article, we turn our attention to correlations of Dirac eigenvalues, i.e. the fluctuations 
about the average behavior of the spectrum. For this purpose we consider multi-level correlation 
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functions. We focus our study on the average connected two-point spectral correlation function, 
Pc(Ai, A 2 ), denned by 



Pc(Ai, A 2 ) = (53«J(Ai - Xi)5(X 2 - A i )) QC D 

~(E 6 ( Xl ~ A fc)>QCD ^( A 2 - Afc))QCD 

k k 

= (5(Ai-A 2 )(^5(A 2 -A fc )) QC D + i?(Ai ! A 2 ). (1.4) 
k 

where (• • -)qcd denotes the average over the QCD partition function, and the are the eigen- 
values of the Dirac operator. In the last line of this equation we have decomposed the correlation 
function into a term containing the self-correlations and the two-point cluster function, R(X\, A 2 ). 
Both terms enter in the disconnected scalar susceptibility which is a more natural object in a 
field theory context. It is defined by 

X = y d mi d m2 lnZ QC D 

= ^(E^+^r 1 E(^+ m 2)~ 1 >QCD (i.5) 

k j 

~y(52( iX k + wi) _1 )qcd (E( iA i + m 2) 1 )qcd 

k j 

= — / dXi dX 2 j— — r. (1.6) 

V J-oo J-oo + mi)(zA 2 + m 2 J 

Because of the averaging over the QCD partition function, p c (Ax,A 2 ) depends on the quark 
masses. Therefore it is not possible to derive the two-point correlation function by inverting 
the relation ( |1.6| ). In order to invert this relation and to compute the correlation function one 
has to introduce special scalar sources, unequivocally related to the eigenvalues of the QCD 
Dirac operator, in the QCD partition function. This partition function contains extra degrees of 
freedom: one fermionic and one bosonic ghost-quark for each special scalar source related to an 
argument of the spectral correlation function. The pqQCD partition function with two special 
sources required for the calculation of the spectral two-point function is given by 



J 1 iu det{i$> + Zl - ji/2) deb(il/> + z 2 - j 2 /2) f} x yv fJ 
Itj(T1i-\ , n\ ^Afc + z\+ ji/2 iX]~ + z 2 + j 2 /2 ' 

= \n ( n . + "»/>) iXt + n _ jl/2 iXt + Z2 _ h/2 

> K J— 1 ' YM 



(1.7) 



Notice that for j% = j 2 = the pqQCD partition function reduces to the QCD partition function. 
Therefore, the derivatives of the pqQCD partition function with respect to the source terms at 
ji = j 2 = are given by averages over the QCD partition function. This enlarged partition 
function was first introduced to study the quenched limit in QCD and is therefore known as 
partially quenched QCD (pqQCD) |Q. It has already been used to compute the spectral density 
of the QCD Dirac operator in j?], g|. In that case, only one such special scalar source had to be 
introduced. 
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We will be interested in the disconnected scalar susceptibility denned by 



X{zi,z 2 ) = — d h d^ InZpqQCD . . 

V Ijl=j2=0 



1.* 



V 



E 



E 



zA/c + 2:1^ zAj + z 2 




QCD 

1 



QCD 



^ iAj + z 2 I 

W J / qcd 



It is related to the spectral two-point correlation function of the QCD Dirac operator (1.4) as 
follows 



1 



+ zi)(iA 2 + z 2 ) 



Because the average over the QCD partition function does not depend on zi and z 2 , this integral 
equation can be inverted. Using that x{ z ii z "i) is °dd in both Z\ and z 2 , we find that 



^Pc(Ai,A 2 ) = ^2 Disc 
1 



zi=i\i, Z2=iM 



X(zi,z 2 ) 



4^2 + e, iA 2 + e) + xM^l + e, iA 2 + e) 

+x(^Ai + e, -iA 2 + e) + x(~^i + e 5 + e) ) ■ 



At low energies, as already discussed in 



(1.10) 



the properties of the spectrum of the QCD Dirac 



operator can be obtained from the low-energy limit of (1.7). As is the case for the usual chiral 
Lagrangian, this effective theory is completely determined by the symmetries of the pqQCD 
partition function. This topic will be discussed in the next section. An alternative way to derive 
the perturbative results for the scalar susceptibility is the use of the replica method [|l4|, 31]. 



However, nonperturbative results cannot be obtained this way [32, 33, 35, 36, 37, [38], 3E] 



Partially quenched Chiral Perturbation Theory, which is the basis of most of our results, 
will be introduced in section 2. The zero momentum sector of this theory will be analyzed in 
section 3. It will be shown that in this domain the two-point spectral correlation function for 
arbitrary topological charge is given by chiral Random Matrix Theory. This nonperturbative 
result generalizes a calculation in [ 40 1 to arbitrary topological charge without starting from chiral 
Random Matrix Theory, but rather from an effective chiral Lagrangian which is obtained from 
the symmetries of the microscopic theory. The contribution of the nonzero momentum modes to 
the scalar susceptibility and its discontinuity is calculated in section 4. In section 5 we evaluate 
the number variance for the different domains in the Dirac spectrum. Concluding remarks are 
made in section 6. 



2 Partially Quenched Chiral Perturbation Theory 

It is well-known that the low-energy limit of QCD is given by a theory of weakly interacting 
Goldstone bosons. The reason is two-fold: the spontaneous breaking of chiral symmetry and 
the existence of a mass gap for the non-Goldstone excitations. In QCD, chiral symmetry is 
maximally broken consistent with the Vafa-Witten theorem, implying that SUi{Nf) x SUn{Nf) 
is broken spontaneously to the diagonal subgroup SUy{Nf). The corresponding low-energy 
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effective theory has been investigated in great detail by means of Chiral Perturbation Theory. 
It describes successfully the strong interaction phenomenology at low energies [41, |4^, 43]. 



In this section we construct the low-energy limit of pqQCD. This low-energy effective theory, 
known as partially quenched Chiral Perturbation Theory (pqChPT), is again solely based on 
the symmetries of pqQCD and is given by a theory of Goldstone modes associated with the 
spontaneous symmetry breaking of chiral symmetry. Because the unitary symmetry is not a 
good symmetry for the bosonic ghost quarks, i.e., unitary transformations violate the complex 
conjugation of the fields necessary to obtain a convergent integral, we start from the complexified 
flavor symmetry group given by Gl(Nf + 2|2) x Gl(Nf + 2|2). The flavor symmetry group of 
the bosonic ghost quarks is then given by Gl{2)/U{2) x Gl(2)/U(2) which results in convergent 
bosonic integrals. The natural choice for the flavor symmetry group of the fermionic quarks 
is U(Nf + 2) x U(Nf + 2). As is the case for QCD we will assume that the axial symmetry 
of the pqQCD is maximally broken by the vacuum expectation value of the chiral condensate. 
We will also assume that supersymmetry is not spontaneously broken. Because of spontaneous 
breaking of chiral symmetry to the diagonal subgroup, the complete Goldstone manifold, which 
also includes fermionic Goldstone modes, is then given by the maximum super-Riemannian 
submanifold of Gl(N f + 2|2) @, @, §. It will be denoted by Gl(N f + 2|2). It consists of a 
fermion-fermion block given by U(Nf + 2) and a boson-boson block given by Gl{2)/U{2). The 
matrix elements of the boson-fermion and the fermion-boson blocks of this manifold are given 
by independent Grassmann variables. The unbroken chiral symmetry group of the pqQCD 
partition function is thus given by Gli(Nf + 2|2) x Gln(Nf + 2|2). 

At low energies the relevant excitations are the Goldstone fields parameterized by 

U = exp(iV2Il a T a /F), (2.11) 

with T a the generators of the Goldstone manifold Gi(N f + 2|2). Under a GIl{Nj + 2|2) x 
Gln(Nf + 2 1 2) transformation of the quark fields, the Goldstone fields transform as 

U -» UlUUr 1 . (2.12) 

The low energy effective partition function is obtained from the requirement that its transforma- 
tion properties are the same as for the pqQCD partition function. Under flavor transformations 
the pqQCD partition function transforms as 

Z pqQGD (£/£.M£/^\0 - iloglSdetf/^ 1 ]) = Z mQCD (M, 6), (2.13) 

where 9 is the vacuum angle, and the quark mass matrix is denoted by M. To lowest order in 
the momenta and quark masses, the effective partition function is thus given by 

Z cS (M,6)= [ ^ dU e~ I d4xC{x) (2.14) 

JUeGl(N f +2\2) 

with the effective chiral Lagrangian given by |30], 0] 

C = ^Str^UdpU- 1 ) - |str(A4t[/ + MU' 1 ) + ^(| - 9) 2 . (2.15) 
4 2 lb 

Here, the quark mass matrix is defined by 

M = diag(zi + ji/2,z 2 + j 2 /2,_m, ■ ■ • ,m,zi - jx/2,z 2 - h/2). (2.16) 

Nf 
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The axial anomaly is included through a mass term for the super-r/ field, <E> = —iFStrlnU, 
while integrating over the full axial symmetry group Gl(Nt + 2|2). This term serves as a 
constraint that projects out the flavor singlet channel. The first two terms in fl2.15|) also appear 
in ChPT to lowest order. However, in the case of pqChPT, there are both fermionic and bosonic 
Goldstone modes. Their masses are given by \/2Sm/i ? , \/2T,zi/F, \/2T,Z2/F, V / E(m + Zx)/F, 
^£(771 + Z2)/F, and \/E(zi + Z 2)/F, depending on their quark content. 

The partially quenched effective partition function was first formulated for the supergroup 
U(Nf + k\k) to study the quenched approximation in QCD |^C|| . This formulation is suitable for 
perturbative calculations. However, an effective partition function based on this group cannot 
be used for nonperturbative calculations of the group integral. For example, the supertrace leads 
to the appearance of both positive and negative masses. The correct integration manifold is the 
super-Riemannian manifold Gl(N f + 2|2) (4|, |, g. 

As already mentioned in the introduction, it is possible to distinguish different domains 
in the spectrum of the Dirac operator. First, the effective partition function is only valid for 
z\ , Z2 "C Aqcd • An important scale is given by the Thouless energy defined as the quark mass 
scale for which the Compton wavelength of the lightest corresponding Goldstone mode is equal 
to the size of the box, i.e. 

m r S f 

■W = J? 

For \zi\ <C m c the Zj-dependence of the condensate is determined by the fluctuations of the 
zero momentum modes. In this domain, the partition function factorizes into a zero momentum 
sector and a nonzero momentum sector [[/], 45, 46]. In |7|, || it was shown that in this domain the 



pqQCD partition function for the one-point function reduces to the chRMT partition function. 
In this article, we will show that the same is true for the spectral two-point correlation function. 
Inside this domain, we can distinguish a second scale: the smallest nonzero eigenvalue 

7T 

Amin = ^y- (2.18) 

For \zi\ S> A m ; n the group integrals can be evaluated perturbatively, whereas for smaller values 
of the Zi the group integrals have to be calculated exactly. In the domain zi > F 2 /SL 2 , the 
nonzero momentum modes become important and have to be taken into account. These three 
domains for the spectral two-point function will be analyzed in detail in the remainder of this 
article. 



3 Quantum and Ergodic domains 

In these domains, corresponding to \z%\ <C m c , the zero momentum mode sector and the nonzero 
momentum sectors of the pqChPT partition function factorize j/J. The Zj-dependence of the 
partition function and therefore the spectral two-point function can be obtained from its zero 
momentum part only. It is given by the super-unitary matrix integral 

ZeS (M,e)=[ „ dUe ^Str(MlU + MU-i)-^±-e)^ (3lg) 

JU€Gl(N f +2\2) 

We decompose the partition function according to 

Z eS (M,9) = [ dve ive Z v {M), (3.20) 
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with the partition function in the sector of topological charge v given by (for an in depth analysis 
of the normalization factors of the different topological sectors we refer to [fl7|]) 

Z V {M) = C ie ~^v [ ^ dUSdet^Ue^^+^^K (3.21) 

JU€Gl(N f +2\2) 

Exact integrals over the Goldstone manifold Gl(N f + 2|2) are mathematically quite complicated. 
In the following we will restrict ourselves to the simplest case, and calculate the scalar suscep- 
tibility in the quantum and ergodic domains only in the quenched limit. A calculation of the 
two-point correlation function based on the supersymmetric formulation of chiral Random Ma- 



trix Theory [48, 49] was first given in [|40|] for the sector of zero topological charge. In this article 
we start directly from the partially quenched chiral Lagrangian and extend the calculation to all 
values of the topological charge. Remarkably, the /c-point correlation functions can be expressed 
in terms of the zero momentum part of the ChPT partition function with 2k additional flavors 
pC(|. The relation of these additional flavors with the equal number of additional flavors that 



enter in the pqQCD partition has been worked out in detail for the one-point function Q. 
3.1 Scalar Susceptibility in the Quenched Limit 

In the quenched limit and in a sector of topological charge the zero-mode part of the effective 
partition function can be written as: 



ZZt(z 1 ,Z2,j 1 ,j 2 )= L DUSdefU exp{— StrM(U + U' 1 )}, (3.22) 

JGl(2\2) 2 

where the mass matrix is given by 

= ( A ? ° ) , (3.23) 



with Mi = T,Vdiag(zi + ji/2,Z{ — ji/2). The Goldstone manifold Gl(2\2) is the maximum 
Riemannian submanifold of G/(2|2) with fermion-fermion block given by U{2) and boson-boson 
block given by Gl(2)/U(2) @, 0, g. 

3.1.1 Parameterization of the Goldstone manifold 

We parameterize the Goldstone manifold in terms of Goldstone modes related with the one-point 
functions corresponding to z\ and , and Goldstone modes that describe two-point correlations. 
A convenient parameterization is given by [40], 



U ~ { w 2 J { -w v 7 !^^ ) { w 2 ) ' [ ' 

where W12 S G7(l|l) and w, w G Gl(l\l). An explicit parameterization of these supermatrices 



will be given below. The advantage of the parameterization ( 3.24 ) becomes clear upon con- 
sideration of the supertrace that appears in the exponent in the effective partition function 

(PD : 



^StrMiU + U- 1 ) =Si +S 2 , (3.25) 
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with 



and 



Si = -Str 
2 



S 2 = -Str 



Mi 



(wiVl — WWW\ + f« 1 — www l 



Miiwi^/^ — WWW2 + w 2 1 \/l~ 



WWWr, 



(3.26) 



(3.27) 



Furthermore, the source terms ji and j'2 only occur in combination with wi and 102, respectively. 
As we will see in the next subsection, the invariant measure of G7(2|2) in the parameterization 
Q3.24 ) factorizes according to 



/J,(w,W,Wi,W2) = fJ l (w,w)fJ,(wi)fJ,(w2)- 



(3.28) 



We thus find that the integrals over w\ and W2 in the computation of the scalar susceptibility 
factorize 

1 8 2 

X(Z1,Z2) 



V dj 1 dj 2 



logZ 



1 



=0 V 



dwdw n(w , w) Ii(w , w)I 2 {w , w) , 



where the integrals Xi(w,w) and l2(w,w) are given by 



dwi fj,(wi) dj^i e" 



3i=0 



1,2. 



(3.29) 



(3.30) 



A further simplification arises by using a polar decomposition of the 2x2 supermatrices that 
appear in the parameterization ( |3.24| ), 

Wi = ViAiv' 1 , i = l,2, (3.31) 
w = vSu^ 1 , w = uSv^ 1 , (3.32) 

where Ai^ and S are 2x2 diagonal supermatrices with commuting elements given by 

A; = diag(e*/ 2 ,e s '/ 2 ), i = l,2, (3.33) 

S = diag(sin#e ip ,isinh</>e i,J ), (3.34) 

S = diag(sin0e- ip ,isinh0e- i<T ), (3.35) 

C = Vl - SS = diag(cos(9,cosh0), (3.36) 

and u, v, vi, and V2, all elements of U(l\l)/U(l) x U(l), can be conveniently parameterized 
according to 



Vi 



exp 



exp 



C 
x 

& 



d = exp 



i = 1,2. 



e 

77 



(3.37) 



After these "Grassmann diagonalizations" , the supertraces in ( 3.26j ) and ( |3.27] ) are given by 



Si = -Str v^vCv^v^Aiv^MiViAi + A^v{ 1 M 1 v 1 A^ 1 



and 



5 2 = -Str 



4 mC« 1 ^ 2 (a 2 u 2 1 M 2 u 2 A2 + A 2 ^2 1 M 2 'U2A2 1 



(3.38) 



(3.39) 
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3.1.2 Measure 



The parameterization of the Goldstone manifold is of the form 

U = WTW = W 2 W~ l TW = W 2 T'. 



(3.40) 



The Berezinian of the transformation from the variables T' to T is one. To calculate the invariant 
measure we thus consider 



T'[U~ 1 dU]T' = W~ 2 dW 2 + \dT'] T 



'1 T-i/ — 1 



(3.41) 



The measure thus factorizes into a product of one factor depending only on W and another factor 
depending only on T" (and on T after the transformation from T" to T). The VF-dependent part 
of the measure trivially factorizes into a wi-dependent piece and a ^-dependent piece. We thus 
find 



d/j,(U) = w l 2 dw\w 2 2 dw\T 1 dT 

= fi(wi)dwi fi{w2)dw2 fi(w,w)dwdw. 



(3.42) 



The first two integration measures simply follow from the invariant measure of G7(l|l)/[/(l|l) 
whereas the integration measure of the T-integrations is given by the invariant measure of 
G/(l|l). Both measures will be calculated in the next part of this section. 

The matrices Wi in the coset G7(l|l)/i7(l|l) have four independent parameters and can be 
parameterized as in (|]H |]37]). We first calculate the measure w, dwi. To obtain the measure 

dw 2 we only have to replace the diagonal elements by their square. The invariant measure 
is given by the Berezinian of the transformation from variables 

to variables dipi, dsi, d^i and dr\i. One easily derives that 

Sw'i = A^Sw" = kr x [8vihi + dAi - AiSvi], 



V; [Wi l dwi]vi 



(3.43) 



(3.44) 

where 5vi = v~ 1 dv{. The Berezinian of the transformation from the variables 5w'( to the dAi 
and the offdiagonal elements of 5i>i is given by 



Sdet 



/ 1 

10 

e s -l 2 - e^/ 2 

\ e*/ 2 - e s ^ 2 

1 



i/2)' 



(3.45) 



The Berezinian of the transformation from 5w" to Sw' is one (factors from the bosonic integra- 
tions cancel against factors from the fermionic integrations), and (5v) 12(^)21 = d^drj. For the 
Berezinian, denoted by B, we thus find 



w i l dwi = Bdipidsid^dr/i 



de^^de^^d&drn 



idsidipid^idrji 



(3.46) 
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The integration measure for the Wi variables is then simply obtained by squaring the diagonal 
elements of Aj , 



H(wi)dwi 



idsidipid^idrji 



e #»)( e - 



(3.47) 



Next we calculate the invariant measure of G7(l|l). This group has eight independent pa- 
rameters and can be parameterized as in ( PT24D and fljT37p . With the observation that ||, |34| 
T~ x dT = dwdw we calculate the integration measure starting from 



dw' = v dwu = SvS + dS — S5u, 

did' = u~ l dwu = 5uS + dS — Sdv, (3.48) 

with 5u = u~ 1 du and Sv = v dv. The Berezinian of the transformation from the variables on 
the left hand side to the variables on the right hand side is given by 

1 1 



(sin 9 + sinh (ft) 2 (cos 2 9 — cosh (ft) 2 

The Jacobian for the transformation of the variables {Su, S22, Sfi, 
is simply given by 

i sinh 2(ft sin 28. 

For the invariant measure of G7(l|l) we thus find 

i sinh 2(f) sin 26 



(3.49) 

S22} to the variables {9, (ft, p, a} 

(3.50) 



T dT = p(w,w)dwdw 



(cos 2 9 — cosh 



-d9d(f>dpdad(dxd£ t dr]. 



(3.51) 



3.1.3 Integration over the Goldstone manifold 

The first supertrace S\ in ( 3.38 ) appearing in the integral ( |3.29| ) can be written as 



Si_ 



( z i + 2-71 ) cos ® cos ^1 ~~ ( Zl 



1 



ji) cosh (ft cosh si 



(3.52) 



+ ji (cos 9 cos ^1 — cosh (ft cosh Si)£i7?i 

+ (cos<9 - cosh0)((zi + -ji) cosV'i - (zi - -Ji) coshsi)(£ - ^i)(r? 

+ ji cos cosh si(cos — cosh (ft) — £1)771 + ^1 (77 — 771)^ 
+ ji(cos# — cosh (ft) (cos ipi — cosh si) £77^1 771. 

The Grassmann integrals in (|3.2S| ) are calculated by collecting the coefficient of £^£1771 i n the 
expansion of e^ 1 . Only the terms linear in j\ contribute to susceptibility. Using the nilpotency 
of the Grassmann variables, such as for example — ^i) 2 = 0, one easily shows that 



— j d£dr}d£,xdr]id h e Sl \ h , 



Q = e EV Z i(cosflcos^-cosh^coshs)( cos0 _ cosh ^( cos ^ i _ coshgl ) 

x (1 + YjV z\ cos 9 cosh ift\ — cosh (ft cosh s\). (3.53) 



An analogous result is obtained from the expansion of e S2 . The integration over the Grassmann 
variables in the disconnected scalar susceptibility ( |3. 29 ) is now trivial. Taking into account the 
measures ( 3.51 , 3.47 ) we arrive at 

»1 



X (zi,z 2 ) = 4E 2 V J dt J 



dp 



tp 



(t 2 -p 2 ) 2 



F(YiVzi)F(£Vz 2 ), 



(3.54) 
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where 



/■27T POO 1 

F[ X ) = (t-p) dtp ds—, - n x(tcosi,-pcoshs) e u(i^-s) 







x (cos ip — cosh s) (l + x(t cos ip — pcosh s)j . (3.55) 

To avoid Efetov-Wegner terms [51, |52| and problems related to the singularity in the integrand 
of T(x), we compute 



1 / 3 h \ t Z 1 ^ 71 " Z" 00 

-^(V + tHW) = — 7 / #/ d s e^^V-Pcoshs) e ^#- S ) 

(i -p) V<9t <9p/ 4 Jo Jo 

x (1 _ e *+#)(i _ e s ~^) [2 + x(t cos i\) - p cosh s) 



(|- + |-) ((* + p)Iv(xt)K u (xp)) . (3.56) 



-<9i dp 

In the first equality we have used the identity 

cos ip — cosh s 1 



(1 _ e -«-*V') j (3.57) 



and the second equality follows from relations for modified Bessel functions. We thus find that 

F{x) = {t 2 - p 2 )I v {xt)K v {xp). (3.58) 

Inserting this expression into ( |3.54| ), the final result for the disconnected scalar susceptibility is 
given by 

X(z u z 2 ) = 4Z 4 V 3 Zl z 2 j dttI„(ZVz 1 t)I v (ZVz 2 t) J dppK u (ZVz lP )K u (ZVz 2 p) 

= ff VZ \? 2 (ziIu+i(^Vz 1 )I v (^Vz 2 ) - Z2l v +i(^Vz 2 )I v {^Vzi)) (3.59) 



x \z 1 K v+1 {Y.Vz l )K v {T.Vz 2 ) - z 2 K v + 1 {Y,Vz 2 )K v {VVz 1 )). 
In the limit Z\ = z 2 , this expression coincides with the result obtained in |l 

3.2 Two-point correlation function 

Finally, the two-point spectral correlation function is given by the discontinuity of the discon- 
nected scalar susceptibility across the imaginary axis CToj) . This follows by using relations for 
Bessel functions such as 

K v {iz) = ^-e m ^ 2 (J I/ (-z)+iN u (-z)), 



K u (-iz) = eT™K v (iz) -ml u (iz), 



TT 



I„(iz) = e-™' 2 J u (-z), -7r<argz< 2 . 
-— = J„(z)N u+1 (z) - J u+1 {z)N v {z). (3.60) 

HZ 



In the quantum and ergodic domains, we then obtain the two-point correlation function 

Pc T ( ^ 2) = *(Ai-A2)y[j£(EVAi)- J, + i(SyA 1 )J i ,_ 1 (SFA 1 )] (3.61) 

AiJ I/+ i(EyAi)j v (syA 2 ) - \ 2 j v +x{iw\ 2 )j v (y.v\{) 



(A? - A|)2 
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The first term represents the contribution due to the selfcorrelations of the eigenvalues (see eq. 
(|1.4|)). It comes from the non-trivial e — ► limit in the discontinuity of the susceptibility ( 1.1C| ). 
This results coincides with the chRMT result |5^| . In the ergodic domain, where SVAi^ 3> 1, 
the two-point spectral correlation function reduces to 

Pc(Ai,A 2 ) = --^ 1 +..., (3.62) 

2tt z (Ai - \ 2 y 

and agrees with the asymptotic result for the two-point correlation function of the Wigner-Dyson 
ensembles. 



4 Ergodic and Diffusive domains 

In these domains, corresponding to A m i n -C \zi\ <C Aqcdj the scalar susceptibility can be com- 



puted perturbatively within pqChPT with Lagrangian for 9 = given by (2.15) 

zp2 y ji#2 n 

C cS = —St^Ud^U- 1 ) - -Str(M(U + U- 1 )) + -^$ 2 + -8^8^ 



(4.63) 



The (Nf + 4) x (Nf + 4) matrix U is parameterized as 

U = exp(iV2IL a T a /F), (4.64) 

with T a the generators of the Goldstone manifold Gl(Nf+2\2), and the singlet field is normalized 

as 

$ = -iFStrlnU. (4.65) 
Its mass is related to the topological susceptibility by M 

» 2(u 2 ) 2mS fl . 

m o = i4 s ~w (4 - 66) 

where m has been introduced to simplify expressions below. For fh = m/Nf this expression 
gives the topological susceptibility for Nf light quarks with mass m, 

<- 2 > - ^- w 

We can distinguish two types of Goldstone modes: those that correspond to the diagonal 
generators and those that correspond to the off-diagonal generators. To one-loop order, the off- 
diagonal Goldstone modes do not mix with the super-singlet field Therefore, their propagator 
is simply given by 

G(p 2 ) = (p 2 + M 2 )-\ (4.68) 

where M is the mass of the corresponding Goldstone mode. The diagonal Goldstone modes, 
on the other hand, do mix with the super-singlet mode. It is still possible to diagonalize the 
quadratic form in the Goldstone fields (see f54} , |30| , 27 ]) . This results in the following propagator 
for diagonal mesons in the sector of fermionic quarks (1 < i,j < Nf + 2), 



. 1 {ap 2 + M 2 ){p 2 + M 2 ) 

' " 1 ' " °V p 2 + M | (p2 + M 2) (p2 + M 2.) ((1 + Nfa)p 2 +M 2 + NfM 2^ > 

(4.69) 



12 



where i denotes the flavor of the quarks of one of the diagonal mesons, and j the flavor of the 
quarks of the other one. The masses in the propagator are given by the Gell-Mann-Oakes- 
Renner relation, M| = TLzijF 2 for i = 1, 2, M 2 = M| = 2£to/F 2 for i = 3, • • • Nf + 2. 

We will consider two limits of this theory: M — > oo (i.e. to — > oo), which is the quenched 
case, and Mq — > oo, which is QCD with Nf light flavors of quarks (and a completely decoupled 
singlet field). In the quenched case, the propagator matrix for the "diagonal mesons" is given 

by 



2 1 {ap 2 + M { 

whereas in the QCD limit we find 



a «^-*WTW-V+>W>+W' fa 1 " i,i - 2, (4 ' 70) 



I (-n 2 -l M 2 ) 



Below we only need the propagators for 1 < i,j < 2 so that both limits can be treated at the 
same time by introducing the propagator 

C (J) 5 1 (^ 2 + 2Zm v /F 2 ) f , 

GiAt ] 'p 2 + mi (p 2 + m 2 ){ V 2 + M 2 n r for 1 ^ t '^ 2 - (4 - 72) 

In the QCD limit we have that rj = 1/Nf and m v = m/Nf, whereas the quenched case is 
obtained by the substitution r/ = a and = to. 

4.1 Disconnected scalar susceptibility 

We compute the disconnected scalar susceptibility for spectral quarks with masses given by Z\ 
and Z2- The mass of all Nf other quarks is taken to be equal to to. To one-loop order the 
contributions represented by the following diagrams have to be taken into account: 





Figure 1: One-loop diagrams in pqChPT which contribute to x{ z i: z 2)- The full lines represent cither 
the "standard" propagator of an off-diagonal meson (4.6?) or the propagator of a diagonal meson (4.72) 
(with a cross). The wiggly lines denote the two different scalar sources. 



(p 2 + M 2 ^ 



with M 2 2 



The propagators in the first diagram in Fig. 1 are given by G(p ) — \p -r ^12) 
S(zi + Z2)/F 2 . In the second diagram of Fig. 1, the lines with a cross denote the propagators 
that mix the diagonal mesons with the super-7/, G\2{p) m ( 4.72 ). Such type of contributions 
also essential in the calculation of the resolvent from the partially quenched chiral Lagrangian 
as well ( [p0|, 0, |27|1). The disconnected scalar susceptibility computed to one-loop order within 
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pqChPT is thus given by 

X (z 1 ,z 2 ) = _^£[(V) 2 + 2G 12 ( F 



2\2 



i 



2\2 



+ 2 



{j]p 2 + 2T,m v /F 2 ) 



VF^l(p 2 + M 2 2 ) 2 (p 2 + M 2 1 ) 2 (p 2 + M 2 2 ) 2 



In the notation of lEoTl with 



G r (M 2 



£(r) 

V 



E(p 2 + M 



2\-r 



the scalar susceptibility can be rewritten as 

2(772:1 - m v ) 2 



X{zi,z 2 ) 



S 2 



G2(M 2 2 ) + V ~ "7> } G a (M?0 + ^-"y G 2 (M^ 



2(772:2 - m f 



(^i-^ 2 ) 2 

F 2 2(772:1 - m ri )(r]z 2 - m n ) 
" S (zi - z 2 ) 3 



(*l - ^) 2 



(4.73) 



(4.74) 



(4.75) 



The functions Gi(M 2 ) and G 2 (M 2 ) for momenta in a finite box were analyzed in detail in [55]. 
They are obviously related to properties of the propagator of a free scalar particle at the origin 
(notice that G 2 (M 2 ) = —d^G^M 2 )). For a box with volume L 4 and momentum cutoff A they 
are given by [42, |55| , 



Gi(M 2 ) 
G 2 (M 2 ) 



1 M 2 
^M 2 log ^ + 5 i (M 2 ,L), 

1 M 2 
"I6^ (1 + 1 ° g -JT)+92(M 2 ,L), 



where 



1 POO 

i6vr2 y ^ 



(4.76) 



(4.77) 



and the sum is over a four dimensional lattice of integers. The functions g r (M 2 ) obviously vanish 
in the thermodynamic limit. For 1/AL 2 -C M <C 1/L they dominate the logarithmic terms in 
the propagators and can be expanded in powers of M resulting in 



9i(M 2 ,L) 
g 2 (M 2 ,L) 



1 



+ O(M /L 2 ) 
+ 0(log(ML)). 



M 2 L 4 
1 

M 4 L 4 



(4.78) 



In the domain 1/AL 2 <C Mjj <C 1/L (z = 1,2) the disconnected scalar susceptibility ( 4.73| ) is 
thus given by 



X(zi,z 2 ) 



irr 



+ 



+ .... 



2 z 2 z 2 V ( Zl + z 2 ) 2 V 
Therefore in the quenched limit, using the relation ( [4.66| ), the result reads, 



2\2 



X{zi,z 2 ) 



2 z 2 z 2 £ 2 V 3 



+ 



>i + z 2 ) 2 V 



+ .... 



(4.79) 



(4.80) 
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and in the QCD limit, with m v = m/Nf, one finds 

■2 



X(zi,z 2 ) 



m 



+ 



1 



+ .... 



(4.81) 



2Nj z\ z\ V ' [z x + z 2 ) 2 V 

In the QCD limit, for spectral quarks with equal masses z\ = z 2 = m, we indeed recover the 
ChPT result that can be easily derived from the results given in [56], 

1} 



N 2 f +2 l 



X 



4Nf m 2 V 



+ 



(4.82) 



In the thermodynamic limit at fixed values of the quark masses, the finite volume corrections 
9i,2(M 2 , L) in ( 4.76| ) can be ignored. This results in the susceptibility 

X 2 



X(zi,z 2 ) 



16vr 2 F 2 

(z 1 -z 2 f 
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In 



Z\ + Z 2 

2/i 



((r}Zi - m^lm^zi + z 2 ) + f]Z\{z\ ~ 3^)]) ln ~ ~ Zl 



z 2 )[r] (z{ + z 2 ) - 2rjm v (zi + z 2 ) + 2m v 



z 2 



(4.83) 



where we have defined the scale [i = A 2 i ?2 /2S (compare to ( 4.76[) ). If the two spectral quarks 
have the same mass z\ = z 2 = z, a somewhat simpler expression is obtained, 



16vr 2 F 4 



^ + ^-(l + 2, 2 )ln^ 

6Z* 6Z V / \i 



(4.84) 



In the limit z — ► 0, the disconnected scalar susceptibility is singular (notice that this expression 
is only valid for z S> 1/VE). This is due to the double pole occurring in the neutral meson 
propagator when m v ^ 0. Such singularities are common in pqChPT. They also appear, for 
instance, in the small mass behavior of the scalar radius of the pion in the quenched approxima- 
tion [57]. In the case m^ = 0, i.e. for {v 2 ) = or in the chiral limit, the algebraic singularities 
in (4.83) cancel and only a logarithmic singularity remains. 

The physically more interesting QCD limit of the disconnected scalar susceptibility, x( z ij z 2), 
is obtained from Q4,83j ) by putting m v = m/Nf and r\ = 1/Nf, 

E X2 



X(zi,z 2 ) 



1 

~16^ 2 V / 2 
2 

N 2 ( Zl -z 2 y 



^ z 1 + z 2 z\ + z 2 + 2m? - 2(zi + z 2 )m 



2(i 



N]( Zl -z 2 ) 2 
| (^zf(zi — 3z 2 ) + 4ziz 2 m — (z\ + z 2 )m 2 ^j In 



£1 
H 



z i 



z 2 



For z\ — ► z 2 
simplifies to 



(4.85) 

the apparent singularities cancel, and in this limit the expression ( |4.85| ) 
T? /m 2 4m 



X( z ) 



+ 



16ir 2 N 2 F 4 \3z 2 3z 



(N] + 2) In 



or, in the chiral limit, 



x( z ) 



T?{N 2 + 2) £ 
16tt 2 7V 2 F 4 ' 



(4.86) 



(4.87) 



Finally, in the case Z\ = z 2 = m, we recover the result derived in [56] within ChPT. 
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4.2 Two-point correlation function 



In this section, we calculate the two-point spectral correlation function from the discontinuities 
of the disconnected scalar susceptibility (see (1.10)) obtained perturbatively within pqChPT in 
the previous section. Therefore, our results for the two-point function are only valid within the 
domain of validity of perturbative pqChPT, 



1/SL 4 < A, < Aqcd i -- 

| Al — A2I S> A mir, ~ 1/SL 4 . 



1,2, 



(4.88) 



The second condition arises from the susceptibilities x(±zAi + e, ^1X2 + e) that contribute to the 



two-point correlation function (see eq. (1.10)). In that case, we have zero-momentum modes 
with mass ~ |Ai — A2I and a perturbative evaluation of the susceptibility is only possible for 
SV|Ai — A2I 3> 1. For eigenvalues inside the domain (|4.88j ) satisfying the additional condition 
Aj <C -F 2 /SL 2 , that is the ergodic domain, only the zero momentum modes contribute to the 
scalar susceptibility which is then given by ( 4.79|) . For the connected two-point correlation 
function we then find 



Pe(Al,A 2 ) = -^2 



1 



+ 



(A1-A2) 2 (A1 + A2) 



+ 



(4.89) 



This result is also obtained from a perturbative expansion of the chRMT result, that is the 



expression (3.61) derived from the zero-momentum sector of the pqChPT partition function in 
Section 3.2. It does not depend on any of the low-energy coupling constants that appear in the 
effective Lagrangian 



^63| ) and is thus the same, independent of the number of flavors and the 
quark masses, provided that chiral symmetry is broken spontaneously 

The perturbative result for the two-point spectral correlation function is obtained from the 
discontinuity of ( 4.750 . At finite volume, in the domain (4.85), it is given by 



Pc(Ai, A 2 ) 



1 



2vr 2 (Ai 
E 2 

"27T2F 4 



A 2 ) 2 2vr 2 (Ai + A 2 ) 2 
" (Ai - A 2 ) 2 



p 1 



+ f^(Al 



A 2 ) 



2^2 



+ 



■pr(Ai + A2) 



(p 4 + ^(Ai + A 2 )2)2- 
rA 2 I 



2\2 



S 2 



S 2 x2 
F? A 2 



(P 



+ ||a 2 ) 2 ( p 4 + w 



^A 2 ) 2 



}. (4.90) 



In the ergodic domain, only the contribution of the first two terms has to be taken into account. 
In the diffusive domain, the thermodynamic limit of ( |4.90| ) is given by 



Pc(Ai, A 2 ) 



E 2 V 
' 64vr 4 F 4 



.(A 2 + A|)( ?? 2 (A 2 + A 2 ) 



16^ 4 



-16vr 



r?Tn T? |Ai||A 2 | 
(|Aa| + |Ai|) : 

VA 2 [Af 



+ 4 



1 



6A 2 A 2 



2m 2 ) 



(A 2 - A 2 ) 2 

,2/\4 , c ^2^2 



(A? - A 2 )3 
4 -~ A? rf\l[\ 



m z v (Xt + 6A 2 A^ + A 4 )ln^| 



3A 4 ]ln^ 2 



6X^X1 



3A 4 ] In M 



(4.91) 



This result can also be obtained directly from the thermodynamic limit of the scalar susceptibility 
in the diffusive domain ( 4.83| ). The two-point correlation function is even in both Ai and A2 
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but is not translational invariant. These properties originate from the pairing of the Dirac 
eigenvalues. For massless quarks (or topological susceptibility equal to zero in the quenched 



case) the expression (|4.9l[) simplifies to 




/vSA.-A 2 ) - - 77T - vrT In - \ R J> + , x2 \ 2 ,M Xi + Xir (4.92) 

Al[A^-6A?Al-3Af]ln^ 

where rj = 1/Nf in the QCD case, and rj = a in the quenched case. Notice that in the limit 
Ai — > A2 the terms proportional to rj 2 are regular. For |Ai — A2I <C Ai f 2 the infrared singular 
part of the correlation function ( 4.92| ) simplifies to 



^,^ = -^lo s ^jM. (4.93) 

This result, derived for m v = 0, does not depend on the parameter rj. It is therefore valid for 
QCD with any number of flavors of massless quarks, even zero. It cannot be obtained from 
chRMT. Beyond the energy scale for which the kinetic term in the chiral Lagrangian cannot be 
neglected, also known as the Thouless energy, pqChPT differs from chRMT. This was observed 
earlier in the analysis of the spectral density in |j [?J and will be discussed in greater detail in 
the next section. 

5 Number Variance 

In the study of disordered systems, a frequently used measure of the spectral correlations is the 
number variance of the eigenvalues. It is defined as the variance of the number of eigenvalues in 
an interval that contains n eigenvalues on average. If the actual number of eigenvalues in each 
such interval for the i-th member of the ensemble is given by rii, the p-th moment of the number 
of eigenvalues is given by the ensemble average (raf). With the number variance denoted by 
S 2 (n) we thus have 

n = (m), 

T?(n) = (nf)-(ntf. (5.94) 

(i) 

If we denote the eigenvalues for the i-th member of the ensemble by A^ the corresponding 
spectral density is given by 

Pi (A) = ]r<5(A-A?). (5.95) 

k 

The number of eigenvalues inside the interval [a, b] is equal to 

m = t 'dXpi(X). (5.96) 

J a 

The average number of eigenvalues inside this interval is thus given by 

n= [ b dX( Pi (X)), (5.97) 

J a 
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and the number variance can be written as 



S 2 (n)=/ / dAidA 2 pc(A 1 ,A 2 ), (5.98) 

J a J a 

with the connected two-point correlation function given by 

MAi, A 2 ) = (tt(Ai)/*i(A 2 )) - (p i (Aa))(p i (A 2 )}. (5.99) 

If the average spectral density is not constant on the interval [a, b] there will be a contribution to 
the number variance due to the variation of the average spectral density. We will always assume 
that this contribution has been eliminated by a procedure called unfolding ||]. 
The two-point correlation function can be decomposed as 

( Pi (\M\')) = <5(A - A') £<«5(A - X k )) + £<*(A - \ k )5(X' - A;)). (5.100) 

k k^l 



The connected two-point correlation function can be written as (see ( |1.4[ )) 

Pc (X, A') = S(X - A')p(A) + R(X, A'), (5.101) 

where the two-point cluster function R(X, A') contains only correlations of different eigenvalues 
and is regular for A — > A'. 

For a finite total number of eigenvalues, we obtain by integration of the connected two-point 
correlation function the sum rule 

/ "dAip c (Ai,A 2 ) = 0. (5.102) 

The contribution due to the self-correlations is canceled by the contribution from the two-point 
cluster function. The coefficient of the linear term in the asymptotic expansion of the number 
variance is given by 

2— / dXp c (— X). (5.103) 



dn T,V Jo T,V ' 

where we have used that p c (Ai,A 2 ) = p c (A 2 ,Ai). This quantity is also known as the spectral 
compressibility. We observe that for large n the coefficient of the linear term vanishes provided 
that the correlation function approaches zero faster than 1/A, 1/A', i.e. when the integration in 
the sum-rule (5.102) is convergent without imposing a cutoff on the total number of eigenvalues. 



If this is the case the correct asymptotic result for the number variance is obtained from the 
asymptotic result of the two-point cluster function provided it is regularized such that the sum- 
rule ( 5.102|) is satisfied. 



In the previous sections, the two-point correlation function was derived under the assumption 
that Ai,A 2 <C Aqcd- On this scale, the variations of the average spectral density J?J can be 
neglected. The average number of eigenvalues in the interval [0, A] is given by 

n= / dXp(X). (5.104) 
Jo 

For m v = the spectral density for A <C Aqcd is given by the Banks-Casher relation 

p(X) = + ■■■. (5.105) 

7T 
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The average number of eigenvalues in the interval [0, A], of eigenvalues is thus given by 



T,VA 



n 



+ ■ 



TT 



(5.106) 



and the number variance can be written as 

Jo 

n + 



£ 2 (n) 



/■7rn/EV 

d\x / d\ 2 p c (Ai, A 2 ) 
J 

irn/Y.V rnn/ZV 

dX 1 / d\ 2 R(X U X 2 ). 
o Jo 



(5.107) 



Perturbative calculations of the spectral two-point function are only valid for |Ai — A 2 | S> 
A m j n and do not include the term due to self-correlations. However, they are included in the 
nonperturbative result ( [3.59 ) valid in the ergodic domain and in the quantum domain. Therefore 
the number variance in the quantum domain and in the ergodic domain is obtained by simply 
inserting the non-perturbative result ( |3.61 ) into ( 5.107| ). We are not aware of any analytic 
expression for these integrals. However both in the case of small n and in the case of large n, 
one can derive concise analytic expressions for the number variance. Only the self-correlations 



in ( 3.61 ) contribute to the small- n limit of the number variance, 



E 2 (n) = n + 0(n 2 ). 



(5.108) 



As discussed above the asymptotic large-n result of the number variance can only be obtained 
from the asymptotic result of the two-point correlation function after it has been regularized 
such that the sum-rule ( |5,102| ) is satisfied. In the ergodic domain this is achieved by 



p c (Ai, A 2 ) = p(X)[5(X - A') + 5(X + A')] - ^ 



1 



+ 



1 



(Ai-A 2 ) 2 + a 2 (A! + A 2 ) 2 + a 



+ 



where a is determined by the sum rule ( 5.102| ), 



(5.109) 



2vrp(A) ' 



(5.110) 



In the ergodic domain, we have that A <C 1/L 2 Aqcd so that p(X) is well approximated by p(0). 
Then p(0) drops out of the expression for the number variance, and at leading order we find the 
familiar logarithmic dependence known from Random Matrix Theory, 



1 



£ (n) = log n + 



(5.111) 



It coincides with the result derived from chRMT. Notice that in the bulk of the spectrum the 
coefficient of log n is a factor of 2 larger. 



At finite volume, for eigenvalues much larger than the smallest eigenvalues but in the domain 
of chiral perturbation theory (i.e. uqqd >t»>l), the number variance can be calculated from 
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n 

Figure 2: The number variance for levels below tiqcd computed from pqChPT for Nf — 3 massless 
quarks. The solid curve represents the result from the zero-mode sector of the partition function, or 
chRMT, and the dashed curve shows the perturbative result. The volume of the box is (20 fm) 4 , the 
Thouless scale is n c w 27.5 and «qcd ~ 10 5 . Notice how late the asymptotic log n-behavior is reached, 
and how early the non-zero momentum corrections are visible. A strikingly similar figure has been 
obtained for a disordered metal in l|58|] . 



the expression for the two-point correlation function that includes the nonzero-momentum modes 
( pOD - We find 



S 2 (n) 



1 

2^2 



1 fc 4 + ( «_y 

4vr 2 



log n + — Tj- 2_\, i0; 



+8vr 2 



A: 1 



rjk 2 _|_ 2Sm I? 
4^ + ~T T ~ 



(5.112) 



where n c = m c /A m i n = F 2 L 2 /ir is the dimensionless Thouless scale. We have used that 
2Trk fl /L with a four-dimensional hypercubic lattice with unit cell of length one. 



In the diffusive domain, the thermodynamic limit of (5.112) is given by 



l + 2rj 2 
16vr 4 



n 



n. 



— log 



n 



n QCD 



(5.113) 



where wqcd = /VA m in ~ ^-QCD/A m in- Therefore, for QCD with Nf flavors of massless quarks 
one finds that in the diffusive domain n c <C n <C ™qcd> the number variance of the Dirac 
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spectrum is given by 



S 2 (n) = -^^f-Vlog^. (5.114) 
V ; 16vr 4 N] \nj 6 n QCD V ; 

The number variance computed from pqChPT for eigenvalues below Aqcd is shown in Fig. 2. 

Finally in the ballistic domain, for n » wqcd the eigenvalues of the Dirac operator in 
QCD are correlated as those for a free Dirac operator. The n 2 log n behavior saturates as n for 
n 3> tjqcd where the interactions become weak because of asymptotic freedom. The spectrum 
thus obeys Poisson statistics, and the number variance is given by 

S 2 (n) = n. (5.115) 



6 Conclusions 

We have analyzed the fluctuations of the eigenvalues of the QCD Dirac operator by means 
of the partially quenched chiral susceptibility. The two-point spectral correlation function is 
obtained from its discontinuity across the imaginary axis. The variance, £ 2 (n), of the number 
of eigenvalues in an interval containing n eigenvalues on average is obtained by integrating the 
spectral two-point correlation function. The generating function for the scalar susceptibility is 
given by the QCD partition function with two additional fermionic quarks and two additional 
bosonic ghost quarks with a mass equal to the spectral mass (also known as valence quark mass) 
that enters in the partially quenched chiral susceptibility. For spectral quark masses well below 
Aqcd the low-energy limit of this theory is completely determined by its global symmetries. 

Based on this partition function we have distinguished three important scales in the Dirac 
spectrum, the smallest nonzero eigenvalue, A m ; n = irfEV, the spectral mass for which the 
Compton wavelength of the corresponding Goldstone boson is equal to the size of the box, 
m c = F 2 /£L 2 , and Aqcd- The analogues of these scales are well-known in mesoscopic physics 
where they separate the quantum domain, the ergodic domain, the diffusive domain and the 
ballistic domain, respectively. 

For spectral quark masses well below m c , in the ergodic domain, the kinetic term in the 
partially quenched chiral Lagrangian decouples from the zero momentum part, and only the 
latter part has to be taken into account. In this domain the partition function only depends on 
the number of flavors, the topological charge and the combination mVT, or zVY*. Nontrivial re- 
sults are obtained by keeping this combination fixed in the thermodynamic limit. A perturbative 
expansion of the zero momentum part of the partition function is possible for spectral quark 
masses that are much larger than A m i n . For smaller spectral quark masses we have calculated ex- 
actly the super-integrals that appear in the partition function for the sector of topological charge 
v. Our results for the scalar susceptibility and the corresponding spectral two-point correlation 
function in this domain are in complete agreement with chiral Random Matrix Theory. This 
result, together with earlier work on the microscopic spectral density, provides strong evidence 
that also all higher order correlation functions that can be derived from the low-energy limit 
of the QCD partition function coincide with the results from chiral Random Matrix Theory. 
The number variance in this domain is given by logn/27r 2 , showing that the fluctuations of the 
eigenvalues are strongly suppressed. 

For spectral masses in the diffusive domain, m c <C z <C Aqcd> when nonzero momentum 
modes have to be taken into account, a perturbative evaluation of the partition is justified. 
In this limit for A <C Aqcd we have calculated the scalar susceptibility to one loop in chiral 



21 



perturbation theory. We have found that in the case of massive quarks (including the quenched 
case) that the scalar susceptibility shows a quadratic infrared divergence. At finite volume this 
divergence is regulated by nonzero momentum modes, and thus only appears for valence quark 
masses beyond the Thouless energy. We expect that this divergence will show up much more 
prominently in lattice QCD simulations than its analogue of quenched chiral logarithms. 

The prediction for the behavior of the two-point correlation function in the diffusive regime is 
Pc(Ai — A2) ~ l/(Ai — A2)^ -4 ^ 2 ^9|. In agreement with this result we have found a logarithmic 
dependence of p c {\\ — A2) on |Ai — A2 1 - The corresponding number variance is given by £ 2 (n) ~ 
n 2 logn. The proportionality constant given by (N 2 + 2) / (167T 2 N 2 F 2 L 2 ) provides us with an 
alternativemethod to determine the pion decay constant. 

The behavior of the number variance is very different in the various domains we have defined. 
Deep in the quantum domain, for small n, the interval contains randomly either zero or one level 
and S 2 (n) = n. When the first level is reached, the number variance almost stops increasing. It 
has a logn behavior up to the Thouless scale n c = 1/F 2 L 2 . Therefore the spectrum is quite stiff 
in this domain. Around the Thouless scale, the number variance increases much faster again. 
Well above n c but below uqcd, the level number corresponding to Aqcdj the number variance 
grows like n 2 logn. Finally, for level number much larger than uqqd, the Dirac spectrum is 
in essence a free particle spectrum because of asymptotic freedom. The levels are therefore 
randomly distributed and the number variance is again given by S 2 (n) = n. This behavior is 
typical for a disordered metallic phase in condensed matter physics. 

In conclusion, we have shown that the fluctuations of the QCD Dirac eigenvalues below 
Aqcd hi the phase of broken chiral symmetry are completely determined by the symmetries of 
the QCD partition function. 
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